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As a result of the quantum, wave-like nature of the physical world, a harmonic 
oscillator can never he completely at rest. Even in the quantum ground state, its 
position will always have fluctuations with variance Ax^p = h/ {2muai), called 
the zero-point motion. This physical theory, which originated over 100 years ago 
to describe the observed behavior of fundamental particles like photons, elec¬ 
trons, and atoms, has also recently been shown to describe the motion of micron- 
scale mechanical systems composed of many atoms (1,2). Although the zero-point 
fluctuations are unavoidable, they can be manipulated. In this work, using mi¬ 
crowave frequency radiation pressure, we both prepare a micron-scale mechan¬ 
ical system in a state near the quantum ground state and then manipulate its 
thermal fluctuations to produce a stationary, quadrature-squeezed state. We de¬ 
duce that the variance of one motional quadrature is 0.80 times the zero-point 
level, or 1 dB of sub-zero-point squeezing. This work is relevant to the quantum 
engineering of states of matter at large^ length scales, the study of decoherence 



of large quantum systems (3), and for the realization of ultra-sensitive sensing of 
force and motion. 


In the quantum ground state, a meehanieal resonator has position fluetuations divided equally 
between its two motional quadratures, Xi and X 2 , whieh are related to the position operator by: 
X = Xi cos Umt + X 2 sin Umt- The ground-state fluetuations minimize the uneertainty relation 
given by the quadratures’ non-zero eommutator: (AXf)(AX|) > | ([Ai,A 2 ]) = {h/2mum) ■ 

Given this uncertainty relation, it is, in principle, possible to squeeze the zero-point noise such that 
fluctuations in one quadrature are reduced below the zero-point level at the expense of increasing 
noise in the orthogonal quadrature. More generally, other non-commuting observable pairs can be 
squeezed, and squeezed states have been created and detected in such varied systems as optical 
(4) and microwave (5) modes, the motion of trapped ions (6), and spin states in an ensemble of 
cold atoms (7). Transient quantum squeezing has also been created and observed in the motion 
of molecular nuclei (3) and of terahertz-frequency phonons in an atomic lattice on picosecond 
timescales (9). In this work, we demonstrate steady-state squeezing of a micron-scale mechanical 
resonator by implementing a reservoir-engineering scheme proposed by Kronwald, Marquardt, and 
Clerk (10). This method is closely related to the approach of Cirac et. al. (11) that was recently 
used to produce quantum squeezed states in the motion of trapped ions. (12) 

A major challenge for quantum squeezing of a radio-frequency mechanical mode is that, even 
at a temperature of 10 mK, the thermal occupation and corresponding position fluctuations are 
far larger than the quantum zero-point fluctuations: Ax^ ~ 100 ■ Ax^^ for a 4 MHz resonator; 
quantum squeezin^can only be accomplished by first overcoming this large thermal contribution. 
In contrast, optical modes are found in the quantum ground state at room temperature. Squeezing 
of mechanical fluctuations was first demonstrated far outside the quantum regime by parametri¬ 
cally modulating the mechanical spring constant (13). Since parametric methods are limited to 3 
dB of steady-state squeezing, the occupation factor of the mechanical mode must be well below 
one phonon to achieve squeezing below the zero-point fluctuations. There are many theoretical 


'in this work, we use “quantum squeezing” to refer to any quadrature-squeezed state with fluctuations below the 
zero-point level in one quadrature. 





proposals for surpassing the 3 dB limit to produce quantum squeezing (14-24), and improvement 
over the 3 dB limit has been realized experimentally with modified parametric techniques (25-27). 
Squeezing below the zero-point fluctuations, however, has yet to be achieved. 

Squeezing via reservoir engineering, as described in (11) and (10), has advantages over other 
methods, as it creates a system in which the mechanics relaxes into a steady-state squeezed state 
without the fast measurements and control necessary for feedback. The scheme consists of apply¬ 
ing two pump tones to an optical or microwave cavity parametrically-coupled to a mechanical res¬ 
onator. The pumps are detuned from the cavity frequency, Uc, by the mechanical frequency, 
and the red-detuned pump has an intracavity occupation, n~, greater than that of the blue-detuned 
pump, (Fig. 1 A). This is a similar set-up to one used for a back-action-evading (BAE) measure¬ 
ment of a single quadrature (28), but with excess red power. In the absence of any blue-detuned 
power, the pump configuration becomes a single, red-detuned drive, as in sideband cooling (29). 
In this limit, the fluctuations of both quadratures are damped and cooled, but the final noise in each 
quadrature is limited by the zero-point noise of the cavity and drive. In the other limit of equal 
red- and blue-detuned power, the pump configuration becomes the balanced drives of BAE. In this 
limit, all back-action noise is added to the X 2 quadrature, leaving the Xi quadrature unperturbed by 
the measurement, and neither quadrature is cooled. In between these limits, both quadratures are 
damped by the excess red power, while the backaction noise added to Xi is less than the zero-point 
noise associated with the total damping. By optimizing the pump ratio between these two limits, 
it is possible in principle to produce arbitrarily-large amounts of sub-zero-point squeezing (i.e., > 
3 dB) if the coupling between the mechanics and the squeezed reservoir sufficiently dominates the 
mechanical dissipation rate (10). 

This interplay between damping and back-action can be seen in the equations for the quadrature 
fluctuations of the system described above. Here, we take the cavity (mechanics) to be coupled 
to a bath with thermal occupation (n^) at a rate k ( 7 ™). Eor optical systems, nf- is usually 
indistinguishable from 0, but for microwave systems, a non-zero is commonly observed at high 
pump powers (30, 31). The cavity and mechanics are coupled with a single-photon coupling rate 



qq = ^Xzp. We also define the enhaneed optomeehanieal eoupling rates for the red and blue 
pumps as G± = go and the effeetive optomeehanieal eoupling rate as ^ — G\. The 

linearized interaetion Hamiltonian for this system is given by 


H = -hd^{G+b^ + G_b) - + H.c, 


( 1 ) 


where d) is the eavity photon ereation operator and W is the meehanieal phonon ereation operator. 
In the good-eavity limit {ujm 3> k), and when k 3> 7 m, the quadrature fluetuations are then given 
by 


(AX2) = x 
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In both eases, the first term is proportional to 2n^ + 1 and has a prefaetor that is less than 1 for 
all ^ > 0. This term represents the damping of both quadratures due to the exeess red power. 
The seeond term is proportional to + 1, and is due to the baek-aetion from the mierowave 
field. We see that the baek-aetion is redueed for (AX^) and inereased for (AX|) relative to 
k) (2’^c^ + l), the amount of baek-aetion we would normally assoeiate with the net damp¬ 
ing rate set by G^- This reduetion for Xi makes it possible to reduee (AX^) below 

It is important to note that the quadrature varianees depend on the thermal oeeupations of both 
the meehanieal and eavity baths, whieh ean both be heated by the applied pump power. These 
thermal oeeupations are also evident in the noise speetrum of photons leaving the eavity. For our 
two-port deviee measured in transmission, the output speetrum derived from the Hamiltonian in 
Eq. [^has the form 


[ce] Sq ~\~ ^out^ 


(7m/2)^ + (w - ^ -f -f 1) 

Tic ^out'^m I . iio ’ 


where f[uS\ = G'^+{'^m/2—iuj) (k/ 2—fee), Kout is the eoupling through the output port of our deviee, 
and Sq is the noise floor of our measurement system. To determine the state of the meehanies, we 
follow an approaeh analogous to that used for meehanieal eooling measurements in the presenee 
of heated baths {32-34)-. we fit the output noise speetrum in the presenee of the squeezing tones 










to extract the bath occupations, and then infer the quadrature variances for a mechanical resonator 
coupled to these baths. 

Our device consists of a mechanical resonator of frequency ojm = 27r x 3.6 MHz coupled 
capacitively to a microwave lumped-element resonator with frequency cUc = 27r x 6.23 GHz (Fig. 
IB). The bare mechanical linewidth is 7 ^ = 27r x 3 Hz at 10 mK, and the cavity linewidth is 
K = 2tt X 450 kHz. The optomechanical coupling rate is = 27i x 36 Hz, and the zero-point 
fluctuations have an amplitude of x^p ~ 2.3 fm. Our measurement set-up is shown in Fig. 1C. We 
can measure either the output noise spectrum of the system (spectral response), or we can sweep 
a small probe tone through the cavity to measure the complex transmission (driven response). 
We first calibrate the normalized spectral power at the output of our system using known thermal 
occupations to find go (Fig. 2A). We then cool our sample temperature to 10 mK, giving an initial 
thermal occupation of ~ 50 quanta. We calibrate the enhanced optomechanical coupling rate, 
G-, to the output power of a red-detuned pump by measuring the total mechanical linewidth, 'ytot = 
'jrn + 4G^/k, from the driven response of the mechanics for linewidths much smaller than k (Fig. 
2B). As we increase the pump power further, the mechanical linewidth becomes comparable to the 
cavity bandwidth, and a strong-coupling model of the system is needed to fit the driven response. 
Fits with this model are in excellent agreement to the low-power calibrations, so we see that the 
device is well behaved and performs linearly over a broad range of pump powers (Fig. 2B). This 
agreement also confirms that the dynamics of our system are well-described by the linearized 
Hamiltonian in Eq. [^in the limit of no blue-detuned power. As we increase the red-detuned pump 
power, the mechanics become more strongly coupled to the cold cavity bath, reaching a minimum 
occupation of 0.22 ± 0.08 quanta (Fig. 2C). The blue-detuned coupling rate, G+, is also calibrated 
to the measured output power using the G_ calibration and a measurement of the output gain at 
0 J+ vs. 0 J-. For further discussion of the system calibrations, see the supplementary online text. 

To squeeze the mechanical motion, a red-detuned pump is applied at Uc — oom and a weaker, 
blue-detuned pump is applied at Uc + Um such that their sidebands overlap at the cavity frequency, 
as described above. We then measure the output spectrum in the lab frame at different pump power 


ratios and fit a spectral model to the background-subtracted output spectrum normalized by the 
measured pump power (Fig. 3A). The spectral model includes “bad-cavity” effects that arise when 

^ 1^, and does not assume that the pumps are aligned at precisely Uc ± Um- In the absence of 
these effects, the model reduces to Eq. Since we determine G_ and G+ from our calibrations, 
and measure n to high precision with a driven response, the only unknown parameters are and 
For parameter estimation, we utilize a Bayesian analysis that systematically incorporates 
the uncertainty from both the nonlinear fitting and the uncertainty in the independent calibration 
measurements. Given the observed data and our knowledge of the system Hamiltonian, we directly 
sample the posterior probability distribution of all system parameters. We then use this distribution 
to extract estimators for the bath occupations, and and for the mechanical quadrature 

occupations, (AXf) and (AX|) (see the supplementary online text for further information). The 
extracted bath occupations are plotted in Fig. 3C, and the quadrature variances are shown in Fig. 
3D. We find that, at our lowest point, we are squeezed to 0.80 ± 0.03 or 1.0 ± 0.2 dB below 
the zero-point fluctuations when n~ = 1.26 x 10^ and = 0.51 x 10^. The spectral fit for this 
point is shown in Fig. 3B. 

To directly measure the fluctuations of a single quadrature, it is possible to introduce a set 
of weak back-action evading tones in addition to the squeezing tones (35) (Fig. 4A). In order 
for the presence of the probe tones not to interfere with the measured mechanical quadrature, the 
probe sideband must be separated from the squeezing sideband by many mechanical linewidths. 
For this device, where the mechanical linewidth at optimal squeezing is comparable to the cavity 
linewidth, such a separation would place the probe sideband outside the cavity bandwidth, sig¬ 
nificantly decreasing the probe gain and leading to unrealistic measurement times. However, we 
have implemented such a measurement scheme using a similar device with a larger cavity band¬ 
width, K = 27r X 860 kHz, and narrower mechanical linewidth jiat = 10 kHz (Fig. 4B). As is 
shown in Fig. 4C, the phase-dependent BAE measurement shows the expected squeezed thermal 
state produced by the two squeezing pumps, with minimum fluctuations of 1.09 ± 0.05 x^p, limited 
in a similar way as the first device by heating. This measurement demonstrates the quadrature 


squeezing produced by imbalanced pumps at ojc ± ujm as described by the theory. 

While a squeezed thermal state always has a positive Wigner function, when the fluctuations 
in one quadrature are reduced below the zero-point level, the squeezed state no longer has a well- 
behaved P-representation - that is, it cannot be represented as an incoherent mixture of coherent 
states (36). For this reason, a quantum squeezed state is considered a non-classical state (37). It is a 
current goal to study the non-classical behavior in larger and larger systems, and recent progress in 
the field of opto- and electromechanics has resulted in the generation of mechanical Fock states (1) 
and entanglement (38); quantum squeezing in a micron-scale mechanical resonator is an important 
addition to this short list. Future work could study the quantum decoherence of a mechanical object 
through the decay of the squeezed state (3). 
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Fig. 1 : Measurement set-up and deviee. (A) Sehematie showing squeezing drive tones detuned 
by ojm from the eavity frequeney. The red-detuned tone has a greater amplitude than the blue- 
detuned tone. (B) Optieal mierograph of the deviee. Gray regions are aluminum, while blue 
regions are the exposed silieon substrate. The eenter square is a parallel-plate eapaeitor with a top 
plate that has a meehanieal degree of freedom. The eapaeitor is surrounded by a spiral induetor, 
and eoupling eapaeitors on the left and right provide input and output eoupling to the deviee. (C) 
Mierowave eireuit diagram for squeezing measurements. The red- and blue-detuned tones are 
filtered at room temperature and attenuated at low temperatures so that they are shot-noise limited 
at the deviee. The deviee is thermally-eonneeted to the mixing plate of a dilution refrigerator. 
Signals are amplified and then measured with a speetrum analyzer or network analyzer. 
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Fig. 2: Calibrations and characterization. (A) We measure the integrated meehanieal sideband 
power (Pm) of a weak red-detuned (red) or blue-detuned (blue) pump and normalize by the mea¬ 
sured pump power (Pthru)- This quantity is proportional to the thermal occupation of the meehanics 
times gQ. The differenee in slopes is due to the differenee in microwave transmission at uj- and 
a;+. (B) Enhaneed coupling rate ealibration for a red-detuned drive. Purple eireles: weak-driving 
regime in whieh the meehanieal suseeptibility is a simple Lorentzian. Blaek line: fit to the weak- 
driving regime damping vs. power. Red eireles: strong-driving regime in whieh the meehanieal 
linewidth beeomes eomparable to the eavity linewidth, and G- must be extraeted using a full 
model. Both the linearity of our deviee at high pump powers and the aeeuraey of the model are 
evident in the adherenee of these strong-driving points to the fit line from the weak-driving regime. 
(C) Charaeterization of sideband eooling in the presenee of a red-detuned drive. Turquoise eir¬ 
eles: eooled meehanieal oeeupation. Yellow squares: eavity thermal oeeupation. The meehanieal 
oeeupation reaehes a minimum of 0.22 ± 0.08 quanta before heating of the cavity bath limits the 
oeeupation. 
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Fig. 3: Squeezing measurements. (A) Background-subtracted output spectra normalized by the 
measured red-detuned pump power at blue/red pump power ratios of 0.3, 0.4, 0.5, 0.6, 0.65, and 
0.7, from red to blue. Darker lines indicate fits from Eq. All spectra are taken at the total power 
of n~ +np = 1.76 X 10^. As the ratio decreases, the damping increases, broadening the mechanical 
linewidth and causing the mechanical noise go from a peak to a dip due to noise squashing. (B) 
Close-up of spectrum and fit at /n~ = 0.4. (C) Values of (yellow circles) and (blue 
squares) obtained from a Bayesian analysis of the spectrum. (D) Left: Calculated quadrature noise 
for (AX^) (red circles) and (AX|) (purple squares). The shaded region indicates sub-zero-point 
squeezing. Right: Close-up of the boxed area showing data with Xi fluctuations below the zero- 
point level. The lowest point has {AX‘f)/xlp = 0.806 ± 0.035. The fit for this point is shown in 
(B). 















Fig. 4: BAE measurement of squeezing. (A) Sehematie of pump and probe tone eonfiguration, 
showing the eavity response and sideband speetra for different probe phases. (B) BAE probe 
sideband speetra for different probe phases with n" = 16 x 10®, rip = 3.2 x 10®, and eaeh 
n^obe _ g gg ^ ]^g6 Quadrature varianee as a funetion of phase, as obtained from the area 
under the sideband lorentzians. The speetra shown in (B) are highlighted in their eorresponding 
eolor. The zero-point varianee is indieated with a dashed black line. 














